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Abstract 

Background: Several studies have reported that the tumor necrosis factor-a (TNF-a) -308G/A polymorphism is associated 
with susceptibility to obstructive sleep apnea-hypopnea syndrome (OSAHS). However, these results are controversial and 
conflicting. 

Objective: To evaluate the association between TNF-a-308G/A and OSAHS risk by meta-analysis. 

Methods: Electronic databases, including PubMed, Embase, China National Knowledge Infrastructure (CNKI), Wanfang, and 
Weipu, were searched to identify relevant studies. Data were extracted from the included studies. A model-free approach 
using odds ratio (OR), generalized odds ratio (ORG) and 95% confidence interval (CI) of the allele contrast to assess the 
association between the -308G/A polymorphism and OSAHS risk. Cumulative and recursive cumulative meta-analyses (CMA) 
were also carried out to investigate the trend and stability of effect sizes as evidence accumulated. 

Results: Seven studies including 1369 OSAHS patients and 1064 controls were identified in this meta-analysis. Significant 
associations were derived from the variants of the allele contrast [(OR, 1.78; 95% CI, 1.45-2.18) or (ORG, 2.01; 95% CI, 1.27- 
3.19). CMA showed a trend of an association. Recursive CMA indicated that more evidence is needed to conclude on the 
status of significance. No significant publication bias was found. 

Conclusions: Our meta-analysis suggested that the TNF-a-308G/A polymorphism contribute to the risk of OSAHS. Further 
studies with larger sample should be performed to confirm our findings. 
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Introduction 

Obstructive sleep apnea-hypopnea syndrome (OSAHS) be- 
comes a major public health problem all over the world recently. 
OSAHS is a respiratory disorder characterized by upper airway 
obstruction during sleep, breathing pauses with oxygen desatura- 
tion, and arousal from sleep. The prevalence of OSAHS is 2% in 
females and 4% in males of middle-aged groups [1]. Studies show 
that OSAHS is an independent risk factor for development of 
cardiovascular, thrombotic and metabolic events [2-3]. Inflam- 
mation plays a central role in this pathophysiological process. 
Tumor necrosis factor (TNF)-a is a potential pro-inflammatory 
cytokine elevated on OSAHS patients and exerted multiple 
physiologic effects on sleep fragmentation [4,5]. Also, chronic 
intermittent hypoxia (CIH) impairs endothelium-dependent vaso- 
dilator mechanism by stimulating NF-KB-mediated TNF-a gen- 
eration [6]. Apparently, understanding the interaction between 



OSAHS and TNF-a is essential for developing therapeutic 
strategies. 

Epidemiological studies confirmed that OSAHS is a chronic 
complex disease with strong genetic components [7,8]. Among 
OSAHS-related candidate genes, TNF-a has been studied 
extensively. TNF-a gene is located on chromosome 6 (6p21), with 
several polymorphisms [e.g., -238 G/A (rs361525), -308 G/A 
(rsl800629), -857G/A (rsl799724) and -1031 T/C (rsl799964)] 
in the promoter region has been identified [9] . Among all the gene 
loci, -308 G/A was mosdy studied [10-16]. 

However, the association between TNF-a-308 G/A polymor- 
phism and OSAHS susceptibility still remained unclear, with 
inconsistent results from different studies. Also, each study has 
limitations in small sample sizes and low statistical power. In order 
to shed some light on this controversy, a meticulous meta-analysis, 
including a cumulative and a recursive cumulative meta-analysis 
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56 studies obtained from initial search 



duplicate studies excluded 



45 studies screened 



12 articles assessed for eligibility 



33 studies were excluded : 
abstracts: n=6 

reviews or conference paper: n=1 0 
non-clinic study: n=7 
irrelevant to OSAHS: n=7 
irrelevant to TNF-a : n=3 



8 articles met all criteria 



4 articles were excluded : irrelevant 
to TNF-a polymorphism 



Figure 1. Flow chart of literature search and study selection. 

doi:10.1371/journal.pone.0106270.g001 

[17], was performed for GAS-related variants in the TNF-a gene 
and OSAHS. 

Materials and Methods 

We performed the present meta-analysis according to the 
guidelines of Preferred Reporting Items for Systematic Reviews 
and Meta-Analysis (PRISMA). 

Publication search strategy 

To identify all eligible studies that address the association 
between the TNF-a-308 G/A polymorphism and OSAHS, we 
systematically searched major electronic literature databases, 
including PubMed, Embase, China National Knowledge Infra- 
structure (CNKI), Wanfang and Weipu. The last search was 
updated to March 2013. The following search terms were used: 
"obstructive sleep apnea hypopnea syndrome" or "OSAHS" or 
"obstructive sleep apnea syndrome" or "OSAS" or "obstructive 
sleep apnea" or "OSA" and "tumor necrosis factor" or "TNF" or 
"tumor necrosis factor-oc" or "TNF-a" in combination with 
"gene" or "polymorphism" or "variants" or "alleles". Addition- 
ally, we manually searched for relevant published studies and 
review articles. No language restrictions were applied. 

Inclusion and exclusion criteria 

Included studies met the following criteria for the present meta- 
analysis: (1) evaluation of the TNF-a-308 G/A polymorphism and 
OSAHS susceptibility, (2) case-control or cohort design, and (3) 
sufficient data to determine the genotype distributions. Studies 
were excluded for the following reasons: (1) absence of control 
subjects, (2) non-clinical studies, (3) reviews, abstracts, or 
conference papers, (4) no reported genotype distribution or allele 
frequency data and participants were not adults (Age ^18 years 
old). We accepted the definition of OSAHS provided by authors in 
each included study. 

Data extraction 

Two investigators (Drs. Zhong and Xiong) independendy 
extracted required data from all the eligible studies. The following 
data were collected from each included study: first author's name 
and publication year, country, ethnicity, genotyping method; age, 
body mass index (BMI), gender distribution, and apnea-hypopnea 
index (AHI) of OSAHS and controls; genotype frequency in cases 
and controls. We also extracted the following additional informa- 
tion: case definition, case and control selection. In addition, it was 
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recorded whether the genotyping in each study was blinded to 
clinical status. If a discrepancy existed in data collecting, it was 
resolved by group discussion. If queries or further study details 
were needed regarding the included studies, the respective authors 
were contacted via electronic mail. 

Quantitative data analysis 

This meta-analysis examined the association between TNF-a 
variant and the risk of OSAHS on the basis of the allele contrast 
(mutant type vs. wild type) [17]. The associations were indicated as 
pooled OR with the corresponding 95% confidence interval (CI). 
In addition, the generalized odds ratio (ORG) was calculated for 
each study [18]. The testing of a genetic association using the 
ORG is powerful even in the absence of Hardy-Weinberg 
equilibrium (HWE) [19]. In synthesizing the studies, the random- 
effects (RE) pooled OR of the allele contrast and the RE ORG 
were used [18]. RE modeling assumes a genuine diversity in the 
results of various studies. Then, the estimation of pooled OR 
incorporates the between-study variance. When there is lack of 
heterogeneity, the RE model coincides with the fixed-effects model 
[17], however, the RE is preferable for statistical inference more 
conservative. Thus, in the present meta-analysis, the RE model 
was used. The heterogeneity between studies was tested using the 
Q_-statistic [20]. If PO_was less than 0.10, then heterogeneity was 
considered statistically significant. Heterogeneity was quantified 
with the I metric, which is independent of the number of studies 
in the meta-analysis [21]. Sensitivity analysis, which examines the 
effect of excluding specific studies, was also performed [17]. The 
distribution of the genotypes in the healthy control group was 
tested to determine whether it is in Hardy-Weinberg equilibrium 
(HWE), an indicator of the quality of a study, Studies with the 
controls not in HWE were subjected to a sensitivity analysis [17]. 
Subgroup analyses were carried out by ethnicity. Furthermore, the 
mode of inheritance was estimated using the degree of dominance 
index (h-index) [22,23]. The h-index is defined as the ratio of the 
natural logarithms of the ORs of the two orthogonal contrasts 
(models): the co-dominant and the additive models. A cumulative 
and recursive cumulative meta-analysis was carried out for each 
polymorphism to evaluate the trend of the risk effect (i.e. OR) of 
the allele contrast in time [17]. In cumulative meta-analysis, 
studies were chronologically ordered by publication year, and then 
the risk effect was obtained at the end of each year, that is, at each 
information step. In recursive cumulative meta-analysis, the 
relative change in pooled OR in each information step was 
calculated. All statistical analyses were performed by Stata (ver. 
11.0, Stata Corporation, College Station, TX) and ORGGASM 
software (http://biomath.med.uth.gr) for the meta-analysis. A P 
value<0.05 was considered statistically significant. 

Results 

Search results 

g § After searching all the electronic databases described above, 56 

~Z § articles were obtained after the initial search. Forty-five studies 

~Z <§ were screened after duplicate studies were excluded. Of these, six 

| "5 were abstracts, 10 were review articles or conference paper, and 7 

ij ° articles were not clinical studies. Additionally, seven of these 

"ro 1 studies were not relevant to TNF-oc, three were not relevant to 

g OSAHS, one study was focused on children, and four were not 

g ^ relevant to the TNF-a polymorphism. Finally, seven articles were 

£ " included in our meta-analysis (Figure 1). 
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Study 

E 



OR (95% CI) 



Weight 



Riha (2005) 
Liu (2006) 
Popko (2008) 
Bhushan (2009) 
Almpanidou (2012) 
Guan (2013) 
Li (2013) 

Overall (l-squared = 18.0%, p = 0.293) 

NOTE: Weights are from random effects analysis 



.135 



o 



1.80(1.20. 2.69) 
3.28(1.45,7.40) 
1.25(0.67, 2.35) 
2.40(1.46, 3.97) 
1.78(1.33, 2.39) 
1.70(1.04, 2.79) 
0.97 (0.45, 2.05) 
1.78(1.45, 2.18) 



19.65 

5.87 

9.36 

13.83 

30.39 

14.12 

6.78 

100.00 



74 



Figure 2. Random-effects (RE) model. Odds ratio (OR) estimates with the corresponding 95% confidence interval (CI) for the allele 
contrast. 

doi:1 0.1 371 /journal.pone.01 06270.g002 



Study characteristics 

Table 1 and table 2 summarize the characteristics of all of the 
included studies. The publication year ranged from 2005 to 2013. 
Half of the articles included were from Asia [11,12,13,14]: three 
from China [12,13,14] and one from India [1 1]. The other studies 
were from western countries: one from the United Kingdom [16], 
one from Poland [15], and one from Greece [10]. Regarding 
ethnicity, four studies concerned patients of Asian origin 
[11,12,13,14], three concerned patients of Caucasian origin 
[10,15,16]. For the determination of the genetic polymorphism, 
validated genotyping methods were used in all studies. The cases 



and controls of included studies were mainly middle-aged and 
slightly overweight, although one study did not report BMI in 
controls [15]. AHI of cases ranged from >21 to 57. 5± 16.3 
events/h; the AHI of controls were all below 5 (Table 1). Of the 7 
studies, one study was not in HWE [11]. All studies provided 
unadjusted estimates of risk effect, and therefore, adjusted pooled 
ORs were not calculated. 

Main results, subgroup analysis, and sensitivity analyses 

For the allele contrast, there were 7 articles evaluating the 
association between TNF polymorphism and the risk of OSAHS. 



Table 3. Comparison of different genetic models for -308G/A polymorphism and OSAHS risk in subgroups. 





Subgroup analysis 


No of studies 


No of cases 


No of controls 


I 2 


PQ 


ORG (95% CI) 


All model-free 


7 


1291 


993 


15.6 


0.31 


1.84(1.47-2.32) 


All in HWE 


6 


1187 


890 


14.6 


0.32 


1.77(1.40-2.25) 


Caucasian 


3 


425 


586 


0.0 


0.52 


1.83(1.42-2.35) 


Asian 


4 


866 


407 


48.2 


0.12 


1.92(1.18-3.13) 


1.86(1.18-2.94)* 


All allele 


7 


1291 


993 


18 


0.29 


1.78(1.45-2.18) 


Caucasian 


3 


425 


586 


0.0 


0.58 


1.71(1.37-2.13) 


Asian 


4 


866 


407 


48.8 


0.12 


1.90(1.24-2.93) 



Abbreviation: HWE, Hardy-Weinberg equilibrium. Note: *mean an adjustment for HWE. 
doi:1 0.1 371 /journal.pone.01 06270.t003 
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Study 




ID 


RR (95% CI) 


a (2005) 


• 1 .45 (1 .06, 1 .98) 


b (2006) 


• > 1.61 (1.20, 2.15) 


c (2008) 


' 1.50 (1.16,1.94) 


d (2009) 


' 1 60 (1 28 1 99) 


e (2012) 


' 1 53(1.30,1.79) 


f (2013) 


' 1.53 (1.32,1.78) 


g (2013) 


' 1.50 (1.29,1.74) 


.466 


2.15 



Figure 3. Cumulative meta-analysis for the allele contrast of 
TNF-a variants and the risk of OSAHS. The pooled odds ratio (OR) 
with the corresponding 95% confidence interval (CI) at the end of each 
year-information step is shown. 
doi:1 0.1 371/journal.pone.01 06270.g003 

The results showed no heterogeneity among studies (I 2 = 18. 0%, 
PQ= 0.293), and the random-effects OR (95% CI), 1.78 (1.45- 
2.18; Figure 2). The sensitivity analysis (exclusion of studies with 
controls not in HWE (Bushan) also showed no heterogeneity 
among studies (r= 13.0%, PQ,= 0.332), and the random-effects 
OR was also significant: OR (95% CI), 1.70(1.37-2.10). An 
adjustment for HWD was performed, and the result did not 
change. In subgroup analysis, namely, no heterogeneity existed 
and significant associations were found in Caucasian OR (95% 
CI), (1.71, 1.37-2.13). In Asian, heterogeneity existed and 
significant associations were found (Table 3). 

In the model-free approach, the results showed no heterogeneity 
among studies (I =16%, PQ_=0.31), and the random-effects 
model was used: ORG (95% CI), 2.01 (1.27-3.19). Similarly, the 
result of the sensitivity analysis (exclusion of studies with controls 
not in HWE (BUSHAN) was ORG (95% CI): 1.77 (1.40-2.25), 
with no heterogeneity (I 2 = 14.6%, PO_= 0.32). An adjustment for 
HWD was performed, and the result did not change. In subgroup 
analysis, significant association was found in both Caucasians 
[ORG (95% CI):1.83 (1.42-2.35); I 2 = 0, PQ=0.52] and Asian, 
[ORG ORG (95% CI):1.92 (1.18-3.13); I 2 = 35%] (Table 3). 

The additive and codominant models also produced significant 
ORs (unadjusted and corrected for deviation from HWE). The 
mode of inheritance, after correcting for deviation from HWE, 
was domiance of allele A (h = 0.50). 

Potential bias 

None of the studies included in the meta-analysis stated that 
genotyping was performed blinded to clinical status. In the 
cumulative meta-analysis, the pooled genetic risk effect remained 
significant in the entire period studied. However, the recursive 
cumulative meta-analysis did not show stability in pooled OR 
change (Figure 3 and Figure 4). Thus, there is a trend of an 
association and more evidence is needed to draw a safe conclusion 
on the significance and magnitude of the effect size. Harbold's test 
indicated that there was no differential magnitude of effect in large 
versus small studies for all settings (P>0.10). 

Discussion 

This meta-analysis examined TNF-oc polymorphism and its 
relationship to susceptibility for OSAHS. The strength of the 
present analysis was based on the accumulation of published data 



2 



1. 5 




2006/2005 2008/2006 2009/2008 2012/2009 2013/2012 

Figure 4. Recursive cumulative meta-analysis for the allele 
contrast. The relative change in the random effects pooled odds ratio 
(OR) in each information step (OR in next year/OR in current year) is 
shown. 

doi:10.1 371/journal.pone.01 06270.g004 

giving greater information to detect significant differences. 
Significant associations were found either for the allele contrast 
or for the model-free approach. In addition, the consistency of 
genetic effects across populations from different ethnicities was 
investigated. Main and subgroup analyses in Caucasian and Asian 
for the allele contrast and the model-free for polymorphisms have 
produced significant results. There was consistency of genetic 
effects across different ethnicities, because they were not differen- 
tiated from the overall effect. 

Currently, meta-analyses of genetic association studies are 
usually based on the analysis of various genetic models. However, 
these models are not independent. In addition, there are no prior 
biological such as the dominant, recessive, additive, and codom- 
inant models; however, these models are not independent 
satisfactions for the choice of a specific genetic model, and the 
interpretation of results is complicated when more than two 
models are significant [23,24]. But in our study, the model-free 
approach was used. And it could provide an integrated way to 
evaluate genetic associations by exploiting all available informa- 
tion in terms of clinical status (case-control) and genotype 
distribution. ORG is a metric that quantifies the probability of 
disease given that a disease subject has a higher or lower 
mutational load than a healthy one. The application of ORG 
might overcome the drawbacks of multiple model testing or 
erroneous model specification and make the interpretation of the 
results easier [18]. The model-free approach (ORG = 2.01) 
indicated that for any two participants, diseased and control, the 
probability of being diseased is 2 folder (relative to the probability 
of being non-diseased) given that the diseased subject has higher 
mutational load than the healthy one, those homozygous for the A 
allele to have Subjects who were homozygous for A allele were 
considered to have the highest mutational load, those homozygous 
for G allele to have the lowest, and those who were heterozygous 
to have an intermediate load [19]. 

In exploring bias, there was a lack of a differential magnitude of 
effect in large versus small studies. After the presence of an 
association was demonstrated, the mode of inheritance was tested 
using the h-index. The presence of significant association, when 
the additive model is significant and the co-dominant model is not 
significant, "non-dominance" or "complete additiveness" is 
assumed and h — 0. When the co-dominant model is significant 
(irrespectively of the significant level of the additive contrast), 
"dominance" is assumed and h>0 or h<0, depending on the risk- 
position of the heterozygote. The cumulative meta-analysis for 
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TNF-oc suggested a consistent significance toward an association as 
evidence from published studies accumulated. The recursive 
cumulative meta-analysis indicated instability in the relative 
change in risk effects, and therefore, more evidence is needed to 
claim or deny the existence of an association between the variants 
and OSAHS. 

We did not identify any published genome-wide association 
study in OSAHS; all studies considered here in were candidate 
gene studies which are known to be susceptible to having limited 
replication validity [25]. The meta-analysis was based on 
unadjusted risk estimates as data for possible effect modifiers that 
influence the estimates of associations (e.g. sex, age, clinical 
manifestations) were not provided by the individual GAS. In 
addition, sampling variability in the GAS may be a possible 
confounding factor on the effect of the genetic variants in OSAHS 
susceptibility. In Bushan's study, controls did not conform to the 
HWE, indicating genotyping errors, population stratification, and 
selection bias in the recruitment of control participants [17]. In 
population terms, lack of HWE implies existence of migration, 
selection, mutation, and absence of random mating; however, 
their impact in candidate-gene studies is limited. An adjustment 
for deviations from HWE was performed to avoid overestimation 
of the genetic risk effects [17]. 

Although significant associations were detected in our study, the 
sample sizes were small, restricting the power of the study. 
Nevertheless, findings of our study could be pooled in a future 
meta-analysis of multiple studies, providing more power to detect 

References 

1. Young T, Palta M, DcmpscyJ, SkatrudJ, Weber S, ct al (1993) The occurrence 
of sleep-disordered breathing among middle-aged adults. N Engl J Med 328: 
1230-1235. 

2. Kono M, Tatsumi K, Saibara T, Nakamura A, Tanabe N, ct al (2007) 
Obstructive sleep apnea syndrome is associated with some components of 
metabolic syndrome. Chest 131: 1387-1392. 

3. Can M, Agikgoz S, Mungan C, Bayraktaroglu T, Kocak E, ct al (2006) Serum 
cardiovascular risk factors in obstructive sleep apnea. Chest 129: 233-237. 

4. Kaushal N, Ramcsh V, Gozal D (2012) TNF-a and temporal changes in sleep 
architecture in mice exposed to sleep fragmentation. PLoS One 7: e45610. 
doi: 1 0. 1 37 1 /journal.ponc.00456 1 0. 

5. Gozal D, Scrpero LD, Kheirandish-Gozal L, Capdcvila OS, Khalyfa A, et al 
(2010) Sleep measures and morning plasma TNF- alpha levels in children with 
sleep-disordered breathing. Sleep 33: 319-325. 

6. Wang B, Yan B, Song D, Ye X, Liu SF (2013) Chronic intermittent hypoxia 
down-regulates endothelial nitric oxide synthase expression by an NF-kB- 
dependent mechanism. Sleep Med 14: 165—171. 

7. Rcdline S, 'fishier PV, Tostcson TD, Williamson J, Kump K, ct al (1995) The 
familial aggregation of obstructive sleep apnea. Am J Respir Crit Care Med 151: 
682-687. 

8. Schwab RJ, Pasirstein M, Kaplan L, Pierson R, Mackley A, ct al (2006) Family 
aggregation of upper airway soft tissue structures in normal subjects and patients 
with sleep apnea. Am J Respir Crit Care Med 173: 453-463. 

9. Ali S, Chopra R, Aggarwal S, Srivastava AK, Kalaiarasan P, ct al (2012) 
Association of variants in BAT1-LTA-TNF-BTNL2 genes within 6p21.3 region 
show graded risk to leprosy in unrelated cohorts of Indian population. Hum 
Genet 131: 703-716. 

10. Almpanidou P, Hadjigcorgiou G, Gourgoulianis K, Papadimitriou A (2012) 
Association of tumor necrosis factor-a gene polymorphism (-308) and obstructive 
sleep apnea-hypopnea syndrome. Hippokratia 16: 217-220. 

1 1 . Bhushan B, Guleria R, Misra A, Luthra K, Vikram NK (2009) TNF-alpha gene 
polymorphism and TNF-alpha levels in obese Asian Indians with obstructive 
sleep apnea. Respir Med 103: 386-392. 

12. GuanJ, Yi H, Mcng L, Su K, Yin S (2013) Study of association of TNF-alpha 
gene polymorphism with OSAHS and metabolic syndrome. Chin Arch 
Otolaryngol Head Neck Surg 20: 8-11. Chinese. 



significant associations [17,26]. The multi-factorial etiology of 
common disorders involving complex epistemic and gene- 
environment interactions reduces the likelihood that a single type 
of studies such as gene-candidate association studies could provide 
conclusive inferences. Our study could not evaluate gene-gene and 
gene-environment interactions due to the small number of patients 
identified. These findings need to be studied in a large cohort of 
OSAHS patients in the future. 

In conclusion, through our meta-analysis, we show a strong 
association between the distribution of genotypes and alleles of the 
single-nucleotide polymorphism TNF-a-308 G/A and OSAHS 
susceptibility. However, additional, well-designed studies using 
modern genetic investigative techniques (GWAS) and larger 
sample sizes are warranted to consolidate these findings. 

Supporting Information 

Checklist SI PRISMA Checklist for systematic review and 

meta-analysis. 

(DOC) 

Author Contributions 

Conceived and designed the experiments: HX MS. Performed the 
experiments: AZ XX. Analyzed the data: AZ HX. Contributed 
reagents/materials/analysis tools: AZ XX HX MS. Wrote the paper: AZ 
HX. 



13. Li TZ, Qian XS, Guo RB, Sun BJ, Liu CT (2013) Relation between TNF-a 
polymorphism and obstructive sleep apnea syndrome accompanying hyperten- 
sion. Acad J Chinese PLA Med Sch. 12: 1265-1268. Chinese. 

14. Liu HG, Guan P, Lin M, Xu YJ, Zhang ZX (2006) The relationship between 
tumor necrosis factor-alpha gene promoter polymorphism and obstructive sleep 
apnea-hypopnea syndrome. Zhonghua Jie He He Hu Xi Za Zhi 29: 596-599. 
Chinese. 

1 5. Popko K, Gorska E, Potapinska O, Wasik M, Stoklosa A, ct al (2008) Frequency 
of distribution of inflammatory cytokines IL-1, IL-6 and TNF-alpha gene 
polymorphism in patients with obstructive sleep apnea. J Physiol Pharmacol 59: 
607-614. 

16. Riha RL, Brandcr P, Vennelle M, McArdle N, Kerr SM, et al (2005) Tumour 
necrosis factor-alpha (-308) gene polymorphism in obstructive sleep apnoca- 
hypopnoea syndrome. Eur Respir J 26: 673—678. 

17. Zintzaras E, Lau J (2008) Synthesis of genetic association studies for 
pertinentgene-disease associations requires appropriate methodological and 
statistical approaches. J Clin Epidemiol 61: 634-645. 

18. Zintzaras E (2010) The generalized odds ratio as a measure of genetic risk effect 
in the analysis and meta-analysis of association studies. Stat Appl Genet Mol Biol 
9: Article 21. 

19. Zintzaras E (2012) The power of generalized odds ratio in assessing association 
in genetic studies with known mode of inheritance. J Appl Stat 39: 2569-2581. 

20. Cochran WG (1954) The combination of estimates from different experiments. 
Biometrics 10: 101-129. 

21. Zintzaras E, Lau J (2008b) Trends in meta-analysis of genetic association 
studies.J Hum Genet 53: 1-9. 

22. Harbord RM, Egger M, Sterne JA (2006) A modified test for small-study effects 
in meta-analyses of controlled trials with binary endpoints. Stat Med 25: 3443- 
3457. 

23. Zintzaras E, Santos M (2011) Estimating the mode of inheritance in genetic 
association studies of qualitative traits based on the degree of dominance index. 
BMC Med Res Mcthodol 11: 171. 

24. Zintzaras E (2008) Variance estimation of allele-bascd odds ratio in the absence 
of Hardy- Weinberg equilibrium. Eur. J. Epidemiol 23: 323-326. 

25. Ioannidis JP, Ntzani EE, Trikalinos TA, Contopoulos-Ioannidis DG (2001) 
Replication validity of genetic association studies. Nat Genet 29: 306-309. 

26. Zintzaras E, Lau J (2008) Trends in meta-analysis of genetic association studies. 
J Hum Genet 53: 176-187. 



PLOS ONE | www.plosone.org 



6 



September 2014 | Volume 9 | Issue 9 | e1 06270 



